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ABSTRACT 



The relatively large number of nearby radio-quiet and thermally emitting isolated neutron stars (INSs) discovered in the ROSAT All- 
Sky Survey, dubbed the "Magnificent Seven", suggests that they belong to a formerly neglected major component of the overall INS 
population. So far, attempts to discover similar INSs beyond the solar vicinity failed to confirm any reliable candidate. The good 
positional accuracy and soft X-ray sensitivity of the EPIC cameras onboard the XMM-Newton satellite allow to efficiently search 
for new thermally emitting INSs. We used the 2XMMp catalogue to select sources with no catalogued candidate counterparts and 
with X-ray spectra similar to those of the Magnificent Seven, but seen at greater distances and thus undergoing higher interstellar 
absorptions. Identifications in more than 170 astronomical catalogues and visual screening allowed to select fewer than 30 good 
INS candidates. In order to rule out alternative identifications, we obtained deep ESO-VLT and SOAR optical imaging for the X-ray 
brightest candidates. We report here on the optical follow-up results of our search and discuss the possible nature of 8 of our candidates. 
A high X-ray-to-optical flux ratio together with a stable flux and soft X-ray spectrum make the brightest source of our sample, 2XMM 
J104608.7-594306, a newly discovered thermally emitting INS. The X-ray source 2XMM J010642.3+005032 has no evident optical 
counterpart and should be further investigated. The remaining X-ray sources are most probably identified with cataclysmic variables 
and active galactic nuclei, as inferred from the colours and flux ratios of their likely optical counterparts. Beyond the finding of 
new thermally emitting INSs, our study aims at constraining the space density of this Galactic population at great distances and at 
determining whether their apparently high density is a local anomaly or not. 

Key words, stars: neutron - X-rays: individual: 2XMM J104608.7-594306 - Catalogs 



1. Introduction 

Around ten years ago, seven X-ray bright thermally emitting and 
radio-quiet isolated neutron stars (INSs) sharing similar prop- 
erties were identified in the ROSAT All-Sky Survey data. This 
group is now commonly ref erred to as the "Magnificent Seven" 
(or M7, for simplicity; see [kaberl 2007, for a review) for the 
reasons that they clearly stand apart from the population of 
rotation-powered radio pulsars. In spite of many searches for 
similar objects in the ROSAT data (e.g. Rutledge et al.l l2003t 
IChieregato et al.ll2~005t lAgtieros et al.ll2006l) . no new candidate 
was identified since the discovery of the last member, RX 
J2143.0+0654 (IZampieri et al.ll2001l) . 

A soft blackbody spectrum undergoing low interstellar ab- 
sorption (kT ~ 40-100eV and Nh ~ few 10 20 ctrT 2 ) is com- 
mon to the seven sources, as well as the absence of a non- 
thermal component extending towards higher energies. We note, 
however, the detection in the optical/UV of a possible non- 
therm al power-law compo nent in the spectrum of RX J0720.4- 
3125 (Ka plan et al.ll2003l) . Very faint optical counterparts with 



blue magnitudes of mg ~ 25-28 were detected for several 
sources, implying high logarithmic X-ray-to-optical flux ratios 
of ~4-5. Relatively to radio pulsars, the neutron star spin peri- 
ods are longer and distributed in a much narrower range, P ~ 3 - 
10 s. Six of the sources show sinusoidal X-ray pulsations with 
pulsed fractions between ~ 1% and 18%. The M7 are believed 
to be nearby (d i= 500 pc), as inferred from the distribution 
of the interstellar medium in the line-of-sight and the equiva- 
lent hydrogen column densities measured in their X-ray spec- 
tra dPosselt et alj|2007h . Furthermore, HST parallaxes are avail- 
able f or the two X-ray brightest members (Ka plan et al.l l2002. 
2007 ) with results largely consistent with those estimated by 
Posselt et al. (2007). The seven sources do not show either per- 
sistent or transie nt radio emission to a r ather sensitive limiting 
flux of ~ lOyuJy ( Kondr atiev et al.ll2008l) and are not associated 
to supernova remnants. Howe ver, unconfirmed claim s of detec- 
tion at long wavelengths exist (Malofeev et al. 2007). 



Proper motion measurements in the optical 
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0633(A), and at the Southern Astrophysical Research Telescope, Cerro 
Pachon, Chile. 



three INSs with bri g ht enough counterp arts, Neuhauser 2001; 
i Kaplan et al l 12002; M otch et al.l 120031; iKaplan et alf 120071; 
iMotch et al J 120051; IZane et all l2006l) a nd in X-rays (for RX 
J1308.6+2127. lMotch et a"i]|2008l 120091) have established these 
objects as cooling, middle-aged (~ 10 5 -10 6 yr) neutron stars, 
probably originated in the nearby OB associations of th e Gould 
Belt. Intense magnetic fields B ~ 10 13 - 10 14 G (see lHaberl 
120071; |vin^Kerkwijk & Kaplan] [20071 for reviews) are inferred 
from the presence of one or more broad absorption lines in the 
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spectra of the sources and fro m the spin-down rates detected in 
three cases by timing studies dKaplan & van Kerkwiikll2005allbl ; 
Ivan Kerkwiik & Ka plan 2008j), where it is assumed magnetic 
dipole braking. For these cases, characteristic ages r c h = P/2P ~ 
1.5-4Myr overpredict by a factor of ~2-10 the kinematic 
travel times inferred b y tracing back the trajectories from their 
likely birth places (e.g. lMotch et alj | 2008f) . The spectral absorp- 
tion lines are commonly interpreted as due to at omic trans itions 
in a partially ionised hydrogen atmosphere (e.g. [Lai 2001) or as 
proton cyclotron features (e.g. lZane et alj|200lh . In spite of the 
high B, their low magnetospheric activity allows for a direct view 
of the surface of the neutron star. Atmosphere models and the 
surface emissivity under extreme ph ysical condition s can thus 
be tested in order to derive radii (e.g. iHo et al.1 12007) and even- 
tually constrain the equation of state of neutron star interior, pro- 
vided that distances are known. Unfortunately, the current lack 
of understanding of the surface composition, magnetic field and 
temperature distributions have limited any definite conclusion so 
far. 

The locus occupied by the M7 in the P-P diagram, some- 
what intermediate between radio pulsars and magnetars, sug- 
gests that they could be linked to other populations of INSs. 
In particular, the similar spin periods and intense magnetic 
fields raised the possibility that some of the M7 could have 
evolved from the young er and more energetic magnetar objects 
dHevl & Kulkarnill9 98). Magnetic field decay provides an addi- 
tional source of heating of the neutron star crust and dramatically 
changes the cooling history of neutron stars with B ^ 10 13 G. 
Aguil era et all d2008l) have shown that the range of observed 
temperatures of the M7 is consistent with them being born as 
magnetars with B ~ 10 14 -10 15 G, provided that their current 
ages are around 10 6 yr (which is in agreement with the kine- 
matic ages). On the other hand, if the neutron star magnetic field 
has not changed substantially over its lifetime, an evolutionary 
link with the high B r adio pulsars (HBPSRs) can be suggested 
(e.g. IZane et al.ll2002l) . Relatively to the M7, these objects are 
generally located at much greater distances and show consider- 
ably higher spin-down energy (E). However, if the HBPSRs are 
a factor of ~ 100 younger, then the difference in E can be ex- 
plained within the standard scenario for puls ar evolution as sum- 
ing the usual magnetic dipole braking model (Kaplan 2008). The 
M7 could then be long-period radio pulsars for which the narrow 
emission beam simply does not sweep over the earth. Finally, the 
only rotat ing radio transient (RRA T) discovered so far at X-ray 
energies dMcLaughlin et alj |2007) also shares strikingly similar 
properties with the M7 - namely the position on the P-P dia- 
gram and the presence of a broad absorption line in its thermal 
spectrum. The source is however hotter (kT ~ 140 eV), fainter 
and expected to be more remote, at ~ 3.6 kpc. 

Considering that within 1 kpc the M7 appear in compara- 
ble numbers as y oung (^ few Myr) radio and y-ray pulsars 
(Popov et al. 2003), they may represent the only identified mem- 
bers of a large, yet undetected, elusive population of radio-quiet 
and thermally emitting INSs in the Galaxy. It should be noted, 
however, that the proximity of OB associations of the Gould Belt 
to the solar vicinity is an important factor that might explain this 
apparent local overdensity of INSs with similar temperatures, 
ages and magnetic field intensities. Since the neutron star cool- 
ing is strongly dependent on mass and, to a lower extent, on mag- 
netic field, a scenario where these stars evolved from a common 
progenitor population of massive stars could be considered. The 
discovery of similar sources at greater distancies is then manda- 
tory in order to make any progress towards understanding their 



Table 1. Results of visual screening of the selected 2XMMp 
sources discriminated by selection region 
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II 


III 


IV 


Known INSs 


16 





1 


1 


INS candidates 


1 


5 


24 


2 


Solar system objects 


2 











Extragalactic X-ray sources 


12 


31 


39 


6 


Bright knots in SNRs 


9 


59 


18 


51 


Diffuse emission in clusters of galaxies 


2 


6 


9 





Spurious out-of-time events / edge 


10 


70 


58 


8 


Spurious in wing of bright object 


2 


15 


18 


6 


Sources in CCD gap 


3 


13 


46 


5 


Other 


2 


3 


8 





Total 


59 


202 


221 


79 



population properties and evolutionary links with other groups 
of Galactic INSs. 

Despite the relatively small field-of-view and low sky cov- 
erage of the XMM-Newton Observatory, its large effective area 
and good positional accuracy at soft X-ray energies make it 
ideal to look for fain t INS candidates beyond the Gould Belt. In 
Pires & Motch (2008) we reported on the preliminary results of 
a programme aimed at identifying new thermally emitting INSs 
in the 2XMMp catalogue. In the present paper we discuss the 
nature of a handful of the X-ray brightest of our INS candi- 
dates, selected from more than 7.2 x 10 4 EPIC pn sources with 
count rates above 0.01 s . The paper is structured as follows: 
in Section [2] we describe the methodology which was applied to 
select the 2XMMp sources and restrict the number of INS can- 
didates; in Section |3] the optical data obtained on a subsample 
of our candidates are described together with their analysis and 
results; Section [4] is devoted to the analysis of the X-ray emis- 
sion of these candidates. Finally, the discussion, conclusions and 
summary are given in Sections [5] and [6] The discovery of the 
new thermally emitting INS 2XMM J104608.7-594306, the X- 
ray brightest among our sample, and the discussion of its nature 
as, in partic ular, a cooling or a n accreting neutron star have been 
presented in lPires et al.ld2009l) . 

2. Selection of candidates 

We searched for new thermally emitting INSs similar to the 
M7 in the pre-release of the XMM-Newton serendipitous source 
catalogue, 2XMMpQ. This version of the catalogue, released in 
July 2006, contains source detections drawn from 2400 individ- 
ual XMM-Newton EPIC observations made between 2000 and 
2006. It contains more than 120,000 unique X-ray sources with 
a median flux of 2.4 x 10~ 14 ergs _1 crrr 2 (in the total energy 
band 0.2- 12keV); around 20% of the catalogued sources have 
fluxes below 1 x 10~ 14 ergs~' cirT 2 . Detection positional accu- 
racy is generally <2" (68% confidence radius) and the total 
sky area covered is ~ 285 deg 2 . Despite the small percentage of 
sky coverage (0.7%), the EPIC instruments provide much better 
positional accuracy at soft energies and faint fluxes relative to 
ROSAT, making the XMM-Newton catalogue a powerful tool to 
look for unidentified INSs. In the near future, the ROSAT succes- 
sor eROSITA^ mission should perform the first imaging all-sky 
survey up to lOkeV with unprecedented spectral and angular 

1 http : //xmmssc-www. star . le . ac . uk/Catalogue/ 
xcat_public_2XMMp .html 

2 http : //www . mpe . mpg . de/pro j ects . html#erosi ta 
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Fig. 1. Left: Boundaries of the 4 different selection regions considered in the HRi XHR2 diagram. Right: Positions of the final 
list of INS candidates in the same diagram (filled circles with error bars). The known /?<9SAr-discovered INSs (stars) occupy the 
less absorbed and cooler portion of the diagram. Dotted lines denote soft absorbed blackbodies of different temperatures (50 eV, 
80 eV, 100 eV, 120 eV, 150 eV and 200 eV, from bottom to top) undergoing hydrogen column densities in range 10 19 to 10 22 cm~ 2 . 
Contours show the hardness ratio distribution of quasars from SDSS DR3 (see text). The candidates selected for follow-up optical 
observations are highlighted by open circles. 



resolution, and can potentially incr ease the number of cooling 
INSs up to one order of magnitude dPosselt et al.ll2008l) . 

In order to select candidates, we only considered detections 
in the EPIC pn camera at a limiting count rate of 0.01 s (0.2- 
12keV), corresponding to a total number of -7.2 x 10 4 X-ray 
sources, and then required these to be well detected and point- 
like. Translating into the more precise catalogue parameters, 
we selected sources with a maximum likelihood of detection 
PN_DET_ML > 8 and required a low probability for extension 
(EP_EXTENT_ML < 4), both properties computed by the cata- 
logue pipeline with the SAS0task emldetect. None of the source 
quality or variability flags provided by the catalogue was used in 
the selection and screening process. 

For the remaining ~4.6x 10 4 entries, we then required as 
main selection criteria (/) no correlations (within 3 <x) with the 
USNO-A2.0 optical catalogue and (ii) a soft energy distribu- 
tion. Applying the first criterion led to the rejection of more than 
2. 1 x 10 4 X-ray sources having a positive match with USNO-A2 
entries. For requirement (ii), instrumental response and auxiliary 
files were used to simulate the count rates and hardness ratios 
(HR) of absorbed blackbodies (in a grid of temperature and col- 
umn density), taking into account the different optical blocking 
filters. Although the effective area of the three (thin, medium and 
thick) filters at energies below 1 keV can differ by a factor of up 
to 3, we verified that their influence in the HR diagrams is minor. 

The HR is defined as 



where C, and C,-+i are source counts into two contiguous of the 
five pre-defined energy bands of the 2XMMp catalogue, 0.2- 

3 http://xmm.esac.esa.int/sas 



0.5 keV, 0.5 - 1 keV, 1 -2 keV, 2-4.5 keV and 4.5 - 12 keV. We 
thus selected the catalogue sources with X-ray energy distribu- 
tion compatible (within 1 <x) with that of template blackbod- 
ies of kT < 200 eV undergoing absorptions in range A^h = 
10 19 - 10 22 cm -2 , by defining their positions in HR diagrams. 

The X-ray column densities of the catalogue sources were 
then verified to be less than or equal to the total Galactic value 
in order to discard intrinsically absorbed objects. Sources with 
very large positional uncertainties (rgo > 4", where r^o is the 
source 90% confidence level error circle) or which were too far 
from the optical axis (at off-axis angles 9 > 11') were also 
discarded. After applying these selections, we ended up with 
~ 1000 sources. Of these, a total additional number of 461 X- 
ray sources with matches in one or more of the optical/IR cata- 
logues USNO-B1.0, SDSS, GSC-2, APM and 2MASS was dis- 
carded as well. We note that our main sources of information on 
the possible astronomical content of the X-ray source error cir- 
cles were the cross-correlations computed by the Astronomical 
Catalogue Data Su bsystem. This mod ule of the XMM-Newton 
reduction pipeline (Watson et al. 2009) provides lists of archival 
entries in over 170 archival catalogues having a position consis- 
tent within 3 cr with that of the 2XMMp source. 

So as to more easily classify the 561 remaining sources 
which fulfilled the above selection criteria, we considered 4 dif- 
ferent regions in the HRi x HR2 diagram, based on temperature 
and column density (see left panel of Fig. [T). We note that, 
although this diagram represents our main means of selecting 
INS candidates, we required soft emission in the HR3 x HR4 
diagram as well so as to avoid X-ray sources with significant 
emission above 2 keV. For every source in these regions, we vi- 
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Fig. 2. Histograms showing the relative frequency of the 
2XMMp sources that fulfilled the selection criteria applied in 
order to find new INS candidates. The four plots show the same 
histograms in the different selection regions (see Fig.[TJ. 



sually checked the X-ray and optical imaged and searched for 
possible identifications in astronomical catalogues and databases 
like NED and Simbad. Source selection and sc reening pro- 
cesse s used the facilities provided by the XCat-DB dMotch et alJ 
I2007I) . a database hosting the 2XMM catalogue and its associ- 
ated pipeline products, including the archival cross-correlations. 

Many of the sources turned out to be false detections in ex- 
tended diffuse emission of mostly supernova remnants (SNRs) 
or due to out-of-time events. The high number of spurious de- 
tections among our selected sources reflects the fact that we did 
not use the source quality flag information in the early steps of 
the selection process. X-ray sources located in the direction of 
nearby galaxies (sometimes with an optical or infrared counter- 
part visible in the finding charts) were classified as "extragalac- 
tic" and thus were no longer regarded as potential INS candi- 
dates. Finally, because the spectral properties of some faint X- 
ray sources are not well determined in short exposures, the last 
step in our procedure was to check if any source with no evident 
classification (and thus a potential INS candidate) was spectrally 
harder in other (longer) XMM-Newton exposures. 

In TableQ]we list the results of the final screening of sources, 
sorted by selection regions. The numbers refer to detections that 



might have entered more than once one or more selection regions 
and not to unique source^- It is worth noting that the known M7 
or, more precisely, the pn observations of the M7 included in 
the compilation of the 2XMMp were among the final list of INS 
candidates, providing a test of the efficiency of our procedure in 
selecting thermally emitting neutron stars. 

Overall, we found 59 very soft and low-absorbed sources 
with X-ray emission consistent with blackbodies of temperatures 
kT < 100 eV and column densities A^h < 1 x 10 21 crrT 2 (region 
I in Fig. [JJ). In this region we retrieved the totality (16) of pn 
observations of the M7 carried out until 2006. Two observations 
of Saturn and a number of bright knots in the extended diffuse 
emission of Galactic SNRs (of the supernova SN 1006 and of 
the relatively old and nearby Cygnus Loop) also are included in 
this region. Interestingly, we noted that there were some cases 
(included as "other" in Table [TJ in which the X-ray source had 
no optical candidate but was nonetheless related to (optically 
bright) high proper motion (HPM) stars. For these cases, we 
checked that the position of the X-ray source was compatible 
with the position of the HPM star at the time of the X-ray obser- 
vation. Region I provided 1 INS candidate. 

More absorbed and slightly hotter sources, consistent with 
blackbodies of kT < 120 eV and N H = (1 -5) x 10 21 cirT 2 (re- 
gion II), totalized a number of 202. This selection region pro- 
vided 4 INS candidates (in addition to the INS candidate that had 
also been selected in region I). Most of the sources consisted of 
bright knots in SNRs - the Galactic remnants of Puppis, RCW 
86 and SN 1006 as well as SNRs in galaxies M33, M31, M51 
and the Small Magellanic Cloud - and of spurious detections 
(mostly of bright objects on the edge of the CCDs or sources 
located in the CCD gaps) and of out-of-time events. 

A total of 221 sources were found in region III {kT = 120- 
200eV and N H = (1 -5) x 10 21 cirT 2 ). Relatively to region II, 
a less important fraction of the total number of selected sources 
consisted of spurious detections in SNRs, which in this region 
were mostly extragalactic. Also relatively to region I and II, a 
larger number of extragalactic X-ray sources were found. Region 
III also provided the largest number of INS candidates (22, in 
addition to 2 repeated entries of selection regions I and II). 
Interestingly, we detected in this region the X-ray emission of 
pulsar PSR J 1722-37 12 and the quiescent neutron stars XMMU 
J164143.8+362758 and CXOU J132619.7-472910 in the glob- 
ular clusters Ml 3 and lo Cen. These two sources were included 
as "other" in Table Q] while the pulsar was included as "known 
INSs". 

Finally, the most absorbed sources, with X-ray emission con- 
sistent with kT < 200 eV and N H = (5 - 10) x 10 21 cm" 2 (region 
IV), corresponded to a number of 79. The detections mainly con- 
sisted of knots in the Galactic SNRs Puppis, RCW 86, RCW 89 
and Tycho as well as in some extragalactic SNRs in the Small 
and Large Magellanic Clouds. This region provided no new INS 
candidate - the objects in Table [TJ corresponding to an INS can- 
didate which already had been selected in region II and to the 
pulsar found in region EH. The results sorted by selection region 
can be visualized in the histograms of Fig. [2] 

In short, out of 7.2 x 10 4 serendipitous EPIC pn sources 
above 0.01 s _1 , fewer than 30 candidates met all the selection 



Finding charts from SDSS or created from digitized optical plates. 



5 In other words, a soft source which was observed several times by 
the XMM-Newton, like one of the M7, can enter the same selection re- 
gion more than once, provided that its spectral properties are compati- 
ble. Similarly, one INS candidate can show HR that are consistent with 
more than one selection region, and in this case it is counted accord- 
ingly. 
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Table 2. Properties of the INS candidates selected for optical follow-up 
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Source identification 


RA 


DEC 


2£ 


H - 


rgq^ 


pn CP/ 


Selection 


Optical 




(J2000) 


(J2000) 


(degrees) 


(cm l ) 


(arcsec) 


/ — 1 \ 

(s ') 


Region 


Data* 1 


zaMM J 1040U8. /-59430O 


1 f\ A C AO "7 

1U 40 Do. / 


-59 43 Uo.l 


a /in 
-O.OU 


1.3d X 1U 


1.33 


U.UoU(4) 


TT TTT J\ J 
11, 111, IV 


*, <> 


2XMM J1210I7.0-464609 


12 10 17.1 


-46 46 11.2 


+ 15.52 


8.70 x 10 20 


2.60 


0.027(6) 


III 


* 


2XMM J010642.3+005032 


01 06 42.4 


+00 50 31.3 


-61.79 


3.18 x 10 20 


3.80 


0.020(5) 


III 


* 


2XMM J043553. 2-102649 


04 35 53.2 


-10 26 50.0 


-34.80 


5.80 x 10 20 


2.70 


0.019(3) 


III 


X 


2XMM J03 1459.9-29 18 16 


03 14 59.9 


-29 18 15.5 


-58.42 


1.32 x 10 20 


2.41 


0.018(4) 


111 


X 


2XMM J214026.1-233222 


21 40 26.2 


-23 32 22.3 


-46.95 


3.51 x 10 20 


1.90 


0.0181(20) 


III 


t 


2XMM J 125904.5-040503 


12 59 04.6 


-04 05 02.3 


+58.73 


1.81 x 10 20 


1.90 


0.0129(21) 


III 


t 


2XMM J 125045 .7-233349 


12 50 45.7 


-23 33 47.7 


+39.31 


7.37 x 10 20 


3.30 


0.0122(21) 


III 


t 



" Galactic latitude. 

Total Galactic extinction jDickev & LockmaJl 990 ) . 
Positional error (90% confidence level). 
d EPIC pn count rate in total energy band (0.2-12 keV). 

e Obtained during different observing periods: * ESO-VLT P79 (pre-imaging), ''SOAR 2007 A, *SOAR 2007B, °SOAR 2008A. 



criteria; these are thus intrinsically soft sources not associated 
to any catalogued optical or infrared object and not likely to be 
spurious. In the right panel of Fig. Q] are shown the positions 
of the final list of unique INS candidates in the hardness ratio 
diagram HRi x HR2. The lowest left part of the diagram is oc- 
cupied by the soft, low-absorbed M7 while our candidates, un- 
dergoing higher photoelectric absorptions, move upwards, along 
with the blackbody lines of hotter temperatures. Contour lines in 
the same diagram sho w the hardness ratio dis tribution of quasars 
from the SDSS DR3 dSchneider et al.ll2005l) which have coun- 
terparts in the 2XMMi catalogue with a maximum likelihood of 
detection greater than 20 (total of 796 sources; the list of cor- 
related sources was extracted using the XCat-DB). It is worth 
noting that their location is right above our selection region III. 
We can see already that several of our INS candidates have hard- 
ness ratios compatible with those of the SDSS quasars. 

We list in Table [2] the equatorial coordinates and count rates 
of a subsample of the candidates, which were object of follow- 
up observations in the optical during the years of 2007 and 2008. 
The choice of sources for the follow-up was made by select- 
ing the X-ray brightest INS candidates which were visible from 
the southern hemisphere. Positions are derived by the catalogue 
pipeline with the SAS task emldetect. The 90% confidence level 

error circle on the position is given by rc,o = 2.15 ^jcr 2 + cr^ yst , 

where <x is the nominal error as given by emldetect. The sys- 
tematic error <x syst on the detection position is provided by the 
2XMMp catalogue and varies with the success or failure of bore- 
sight correction of the EPIC sources with optical matches in the 
USNO-B 1.0 catalogue. 

3. Optical follow-up 

We are presently conducting an optical campaign with the pur- 
pose of investigating the nature of the X-ray brightest INS candi- 
dates of our sample. The immediate goal is to find a possible al- 
ternative identification - essentially faint polar-type cataclysmic 
variables (CVs), late-type stars or active galactic nuclei (AGN) 
- to the selected sources, through the analysis of the spectra and 
the colour indexes of the optical objects that may be present in 
the X-ray error circles. 

The 8 INS candidates in Table [2] were observed using the 
8.2 m ESO Very Large Telescope (ESO-VLT) and the 4.1m 
Southern Astrophysical Research Telescope (SOAR) facilities in 
Chile during 2007 and 2008. Table [3] shows the log of all the 



observations that we have obtained and analysed so far. We note 
that the exposure times, seeing and airmasses listed in Table[3]are 
averaged per filter. Hereafter, we adopt the convention "XMM 
Jhhmm" to designate the 8 X-ray sources discussed in this work. 

3.1. Observations and data reduction 

We obtained deep imaging of the fields of our candidates un- 
der photometric sky conditions (see the log of observations 
in Table [3]>. The FOcal Reducer low/dispersion Spectrograph 
(FOR S2:lAppenzeller et all Il998l) and SOAR Optical Imager 
(SOL lSchwarz et al.l2004l) were used. The FORS2 detector con- 
sists of a mosaic of two 2 k x 4 k MIT CCDs (15 fim) and is op- 
timized for the red band. It provides imaging at a pixel scale 
of 0.25"pixeL' using standard read-out mode and a 2 x 2 bin- 
ning (6.8' x 6.8' field-of-view). With this configuration the gain 
is 1.25 e" ADU 1 and the read-out noise 2.7 e". The SOI instru- 
ment uses a mosaic of two E2V 2kx4k CCDs (15 /mi) and it 
is optimized for the blue and UV bands. It has a field-of-view of 
5.3' X 5.3' and a pixel scale of 0.15"pixel _1 using a 2x2 bin- 
ning. It was used in slow read-out mode, minimizing the read-out 
noise to 3.1 e" with a gain of 0.4 e" ADU -1 . 

Only pre-imaging dat£0 were obtained with the ESO-VLT 
while ~ 90% of the proposed observing time was executed by 
SOAR in three different observing periods. Dithering patterns 
with offsets of 2" and 5" were chosen for the SOAR and the 
ESO-VLT observations, respectively. Seeing during the several 
nights of observation varied between 0.58" and 1.2"FWHM, 
with mean values of 0.70" and 0.88" for the ESO-VLT and 
SOAR nights. 

We used the already bias and flat field corrected frames pro- 
vided by the ESO-VLT pipeline in the analysis. Images in each 
filter were combined using an algorithm to remove cosmic rays 
and bad pixels based on CCD statistics (cedclip in IRAff] task 
imcombine). SOAR data were reduce d using standard proce- 
dures within IRAF V2.14 dTodvl [l986). The observations were 
bias and overscan subtracted and combined skyflat and dome 
flat images in each filter were used to correct the pixel-to-pixel 
response variations across the CCDs. The science frames in each 
filter were also combined to eliminate cosmic-rays and increase 
the signal-to-noise ratio of the data. 

6 The pre-imaging had the purpose of selecting potential optical can- 
didates for spectroscopy or deep imaging. 

7 http://iraf.noao.edu 
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Table 3. Description of the optical observations 



Target 


Night 


Telescope 


Exposures 


Exp. Time 


Total Exp. Time 


FWHM 




Airmass 










(s) 


(s) 


(arcsec) 






XMM J 1046 


22-02-2007 


ESO-VLT 


2xB, 2xR 


150, 150 


300, 300 


0.81,0.68 




1.27, 1.26 


XMM J1210 


25-02-2007 


ESO-VLT 


2xB, 2xR 


150, 150 


300, 300 


0.83, 0.70 




1.16, 1.15 


XMM J1250 


09-05-2007 


SOAR 


14 xU, 7xB, 3xR 


500, 200, 200 


7000, 1400, 600 


0.90, 0.74, 


0.80 


1.10, 1.03, 1.01 


XMM J1259 


10-07-2007 


SOAR 


4xU, 4xB, 4xR 


500, 200, 200 


2000, 800, 800 


0.75, 0.68, 


0.62 


1.02, 1.01, 1.01 


XMM J2140 


11-07-2007 


SOAR 


4xU, 4xB, 3xR 


500, 200, 200 


2000, 800, 600 


1.11,0.95, 0.80 


1.19, 1.23, 1.25 


XMM JO 106 


21-07-2007 


ESO-VLT 


2xB, 2xR 


150, 150 


300, 300 


0.62, 0.58 




1.16, 1.15 


XMM J0314 


29-01-2008 


SOAR 


4xU, 5xB, 4xR 


500, 200, 200 


2000, 1000, 800 


0.85, 0.82, 


0.72 


1.18, 1.22, 1.24 


XMM J0435 


30-01-2008 


SOAR 


4xU, 4xB, 4xR 


500, 200, 200 


2000, 800, 800 


1.14, 1.16, 


1.10 


1.10, 1.11, 1.12 


XMM J 1046 


08-02-2008 


SOAR 


6x V 


1120 


6720 


0.94 




1.23 



The observation of photometric standard stars fields dLandoltl 

fl992h was arranged to be executed during the same night as our 
SOAR targets. Fields RU 149, Mark A and PG1323-085 were 
selected to calibrate the observations of sources XMM J0435, 
XMM J0314, XMM J2140, XMM J1259 and XMM J1250. The 
deeper observations of XMM J 1046 were calibrated with the 
standards PG1047+003 and SA 104. For every field (as appro- 
priate depending on the requested observations of the science 
target) we obtained 10 s (U or V) and 5 s (B or R) single expo- 
sures. The airmasses at which the observations of standard stars 
were carried out were similar to the ones of the science targets. 
As only one exposure was taken in each filter (at a given airmass) 
it was not possible to determine the extinction coefficients; mean 
values for the periods were adopted instead. On the other hand, 
the zero-point magnitudes and colour coefficients were deter- 
mined from our observations of the standard stars in each night. 
The fit to the standard system was performed relative to the U-B 
(for the U magnitudes) and B - R (for the B and R magnitudes) 
colours. The V SOAR instrumental magnitudes (period 2008A) 
were transformed adopting a null colour coefficient and the mean 
value of the extinction coefficient for the period. The instrumen- 
tal magnitudes of the ESO-VLT images were transformed using 
directly the extinction coefficients and colour terms provided by 
the ESO-VLT calibration web pagefl We compared the R mag- 
nitudes of non-saturated GSC-2 stars present in the fields with 
our results and concluded that the agreement is good with a mean 
dispersion of ~ 0.03 mag. 

3.2. Data analysis and results 

3.2.1. Astrometry 

The astrometric calibration was performed using the USNO- 
B1.0, 2MASS and GSC-2 catalogues and the GAIA 4.2-1 soft- 
ware£J The pixel coordinates of the non-saturated catalogued 
stars in every combined science frame were determined by fit- 
ting a two-dimensional Gaussian function to their intensity pro- 
files. An astrometrical solution was then computed by fitting the 
pixel to the celestial coordinates. Overall, our astrometric errors 
are of ~ 0. 15" or better. 

3.2.2. Photometry 

Two different approaches were adopted to measure magni tudes: 
stand ard PSF fitting, as implemented in the DAOPHOT (IStetsonl 
|1987|) package for the IRAF environment and, for the fields 



of candidates located at high Galactic latitudes (\b\ ^ 20°, see 
Table |2J, variable elliptic apertures and background maps using 

SExtractor 2.5.C0 

PSF fitting using a two-dimensional Gaussian function was 
adopted whenever there were a fair number of isolated and bright 
non-saturated stars available, allowing a good determination of 
the modeled PSF to apply to all detected sources on the image 
(including faint and overlapping stars). We note, however, that 
none of the fields is severely crowded, not even the one located 
at b ~ 0°. When this was not possible, we adopted SExtractor 
instead. This software is most commonly used to reduce galaxy 
survey data but it also performs well in moderately crowded stel- 
lar fields. As many extragalactic objects were likely to be present 
in our high Galactic latitude fields, its usage is convenient since 
it is possible to derive, for instance, information on how elon- 
gated an object is; additionally, making use of a neural network, 
SExtrator classifies any given object as stellar or non-stellar. 
However, this classification is less reliable for faint fluxes. Its 
main advantage is that flux measurements of extended objects 
tend to be more accurate when compared to the ones obtained 
using fixed aperture photometry, since elliptical apertures with 
variable sizes (based on the object intensity) are used instead. 
Another advantage is that a global background map is created, 
which better accounts for local spatial variations of brightness 
due e.g. to nebular emission or caused by scattered light from 
bright objects. 

Using SExtractor, we created catalogues gathering position, 
intensity and shape information on every optical object present 
in our high Galactic latitude frames. Source detection was car- 
ried out by convolving a Gaussian filter having the mean image 
FWHM. The resulting list of sources in each field was then cor- 
related to keep the optical objects which were detected in all 
filters. In the following subsections we use the information ex- 
tracted from these optical catalogues to discuss the nature of our 
sample of INS candidates and field objects. 

3.2.3. X-ray/optical associations 

We found at least one optical object inside the X-ray error circles 
(90% confidence level) of all INS candidates but one, the X-ray 
brightest source XMM J1046 (see Fig. [3). In this case no optical 
object brighter than the limiting magnitude of our present data 
is present within ~ 4.3" (^ 5 cr) from the position of the X-ray 
source (see discussion in 13. 2. 5b . 

Table |] lists the R magnitudes and B - R and U - 43 
colours of the possible optical counterparts of the remaining 



8 http : //www . eso . org/observing/dfo/quality 

9 http : // star-www . dur . ac . uk/~pdraper/gaia/gaia . html 



http : //terapix . iap . fr/rubrique . php?id_rubrique=91 
11 Only for SOAR data. 
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Fig. 3. R images of the fields of the observed INS candidates. Candidates XMMJ1046, XMM J1210 and XMM J0106 were observed 
with the ESO-VLT in period P79 while candidates XMM J0435, XMM J0314, XMM J2140, XMM J1259 and XMM J1250 were 
observed using the SOAR telescope. In these images, north and east point upwards and to the left, respectively. Circles show the 
90% confidence level error on the positions. 

Table 4. Results of follow-up investigations 



Candidate 


r 


R 


U -B 


B-R 


LR 


^chance 


e 






(arcsec) 










(%) 






XMM J1210 


1.55 


20.13(3) 


7 


1.35(3) 


8.4 


2.8 


0.037 


0.18 


XMM JO 106 


2.75 


24.51(9) 


? 


£2.0 


1.2 


22.4 


0.178 


1.93 




1.68 


24.58(9) 


? 


>1.9 


1.2 


22.4 


0.162 


1.96 


XMM J0435 


1.20 


22.15(4) 


-0.51(8) 


0.93(7) 


15.2 


2.3 


0.068 


0.77 


XMM J0314 


0.67 


21.37(3) 


-0.13(7) 


1.43(4) 


32.1 


1.7 


0.022 


0.41 


XMM J2 140 


1.26 


23.78(9) 


-1.08(16) 


1.06(15) 


6.0 


6.4 


0.156 


1.40 


XMM J 1259 


1.34 


21.64(4) 


-1.08(4) 


0.76(5) 


25.5 


1.0 


0.098 


0.61 


XMM J 1250 


1.80 


22.19(4) 


-0.71(4) 


0.74(4) 


8.4 


4.1 


0.029 


0.57 



The X-ray flux is computed in the 0.15-3 keV energy band assuming an absorbed power-law model. 



seven sources, as well as their positional offsets relative to the 
X-ray coordinates and the implied logarithmic X-ray-to-optical 
flux ratios, uncorrected for absorption. The candidate counter- 
parts were well detected in all (UBR or BR) filters, with the ex- 
ception of the two optical objects possibly associated with XMM 
J0106, only detected on the R image. Among the possible coun- 
terparts, these were also the faintest ones (R ~ 24.5). 

To estimate the reliability of each X-ray/optical association, 
we defined a likelihood ratio (LR) based on the probabilities that, 
given the positions, an optical object is the true counterpart of 
the X-ray source or is just a random o bject that happens to li e 
angularly close to its position. Followin dde Ruiter et alj(ll977l) . 
LR can be defined as: 



iS = s exp 



t("-0 



(2) 



where A is defined as A — ncp-p, cr is the error on the X-ray 
and optical positions and p is the surface density of brighter or 



equally bright objects on the image. The dimensionless variable 
x = r/cr is the significance of the position offset (r). 

We list in Table H] the LR quantities for every possible X- 
ray/optical association, which ranges from ~ 1 to 32. Although 
objects with large LR are more likely to constitute a true asso- 
ciation, the computation of actual probabilities of identification 
requires a calibration of LR taking into account the a priori prob- 
ability that any X-ray so urce has a counterp art in the optical sam- 
ple considered (see e.g. lPineau et al.ll2008l) . 

LR strongly depends on the local density of optical objects; 
the chance of a spurious association rising with the number of 
objects of similar or brighter magnitudes. In order to derive the 
LR distribution of spurious associations and therefore estimate 
the chance probability of X-ray/optical associations in our data, 
we performed Monte Carlo simulations randomly changing the 
position of the X-ray source accross the optical fields. For a 
given simulated X-ray source, we then computed the LR quan- 
tity given by Equation (O for every optical object (if any) lying 
inside ajq. The overall distribution of positive matches (among 
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Fig. 4. Histograms showing the LR distribution of X-ray/optical associations in the fields of candidates for which possible optical 
counterparts were found. Each histogram is computed from an original number of 1000 Monte Carlo simulations. The number of 
matchings is given in the right upper corners. Hashed bars show simulated X-ray/optical associations which have a LR greater than 
or equal to the one found for the X-ray source. 



1000 simulations) can be seen in the histograms of Fig. [4] As 
expected, the distribution of LR is highly asymmetrical, strongly 
peaking at low values and then monotonically extending towards 
high LR. Hashed bars in Fig.|4]highlight the number of simulated 
random X-ray/optical associations with LRj > LRx, where LR, 
and LRx are the likelihood ratios of simulated and the actual X- 
ray/optical association of the given field. 

In general, the simulations show that the probability of a 
chance association is low for the optical objects found inside r^Q 
on the real data (Table 3). Out of 1000 simulations, we usually 
found that less than 3% of the cases had an X-ray /optical asso- 
ciation with a likelihood ratio greater than or equal to the real 
one. The two exceptions are the associations found for source 
XMM J0106, with Prance ~ 22% and, at a smaller degree, for 
source XMM J2140, with f ch ance -6%. These are the two X- 
ray sources with the faintest optical candidates and the only 
ones with log(/x//R) > 1. However, whereas the two candi- 
date counterparts for XMM JO 106 were detected only in the R 
band, the one for XMM J2140 has been well detected in the U, 
B and R filters, and it shows very blue colours. We argue in the 
following that, based on the results obtained here and together 
with the analysis of their colours and X-ray-to-optical flux ratios 
(Section |3.2.4| ), the optical objects found inside the X-ray error 
circles of all INS candidates are very likely to constitute their 
true optical counterparts, with, again, the exception of source 
XMMJ0106. 

3.2.4. Identification with other classes of X-ray emitters 

As mentioned before, some populations of astrophysical objects 
are expected to pollute our sample of INS candidates. In order 
to identify these among the sources having optical candidates, 
we compared their positions in the colour-colour g' — r 1 Xu' - V 



diagram with those of the major classes of X-ray emitters (left 
pannel of Fig. [5]). The magnitudes in the Sloan photometric sys- 
tem were computed usi ng the Johnson-Mo rgan-Cousins trans- 
formation equations in Fukugita et al. dl996l) . In this diagram 
are plotted, as c ontours, the spectro scopically identified popula- 
tion of quasars jRichards et alj|200ll) and cataclysmic variables 
dDillon et al.l 120081 an d references ther ein) from the SDSS, as 
well as late-type stars dBilir et al.| [2005). Only the five INS can- 
didates observed with SOAR are shown, since we do not have U 
magnitudes for the two candidates observed with the ESO-VLT 
(sources XMM J1210 and XMM J0106). We note, however, that 
the B-R colours of the optical candidates of these sources stand 
among the "reddest" of our sample (Tabled. 

The u' - g' and g' - r' colours of the optical candidates of 
the SOAR targets are clearly consistent with those of AGN and 
CVs. Among the five sources, the proposed counterpart of XMM 
J03 14 is the only one which stands somewhat more isolated from 
the bulk of these two populations of X- ray emitters, al t hough 
some quasars at still low redshifts from iRichards et aLl ([2001) 
do show similar colours. None of the optical sources exhibits 
colours typical of late-type stars. The strong blue/UV excess, 
with u' — g' < 0, shown by the proposed counterparts of sources 
XMM J1259 and XMM J2140 is remarkable. The colours of the 
optical objects in the error circles of XMM J 1250 and XMM 
J0435 are somewhat more intermediate but still consistent with 
those of AGN and CVs. 

Interestingly, when compared to other uncatalogued field ob- 
jects present in our SOAR frames and detected in all filters, the 
colours and even the UV excess exhibited by our sample of opti- 
cal candidates do not stand out from the overall field population 
(see right pannel of Fig.|5|. For further comparison, we also plot 
as density contours SDSS objects present in five regions of the 
sky located at similar column densities as the SOAR fields. More 
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Fig. 5. Left: Colour-colour g' —r 1 Xu' - g' diagram comparing the locations of different astrophysical populations of X-ray emitters. 
Red stars (and contours) show the position of identified QSOs at redshifts lower than 3. Hashed blue and solid green contours are 
CVs and late-type stars from the SDSS, respectively (see text). The five INS candidates observed with SOAR are represented by 
plain black circles with labels. Right: Same diagram showing the colours of field objects detected in the three (U, B and R) filters 
of our SOAR observations (plain red circles). The population of SDSS objects located at a similar column densities is shown for 
comparison as a density contour. Open triangles are identified quasars located at redshifts greater than 3. Optical candidates of the 
soft X-ray sources are shown as large blue circles. 



Table 5. Upper limits on the B, V and R magnitudes for the op- 
tical counterpart of XMM J 1046 



Filter 




C dct 


m AA 


log(F x // r . i ) 


B 


26.0 


2.4 


26.2(5) 


£2.2 


V 


25.5 


2.5 


25.9(4) 


£3.1 


R 


24.2 


2.4 


25.5(5) 


£3.2 



precisely, we gathered SDSS entries in a SOAR field-of-view 
centered at five different directions, requiring r' magnitudes be- 
tween 21 and 24, which is the range covered by our SOAR field 
objects, and u', g' and r' errors smaller than 0.5. Overall, the 
plot shows that the SOAR and SDSS objects constitute the same 
population of "blue field objects" located at high Galactic lati- 
tudes, with colours not particularly distinguishable from those of 
the sample of soft X-ray counterparts. Apart from these blue ob- 
jects, a smaller fraction of the SOAR field population occupies 
the same location in the colour diagram as the one of quasars 
identified at more remote redshifts (z > 3; see open triangles 
in Fig. [5). Most (k 99%) of the SOAR objects do not have a 
correlation with the 2X MMi and none with the FIRST (radio) 
catalogues. According to Ma teos et al.l d2008l) . the expected den- 
sity of soft (0.5 -2 keV) X-ray sources at high Galactic latitudes 
\b\ > 20°, with fluxes greater than fx > 10~ 14 ergs~' cirT 2 is 
~ 100 deg~ 2 , corresponding to ~ 0.75 sources in a typical SOAR 
field-of-view - which is roughly in agreement with our results. 
In summary, despite the blue colours of the field SOAR objects 
are similar to the ones of the INS candidates, they are not active 
in X-rays or radio wavelengths. 



We analysed the extent distribution of field objects and opti- 
cal counterparts using the ellipticity derived by SExtractor. The 
goal was to verify the possibility of separating stars/ AGN from 
galaxies by measuring their point spread functions. Defining el- 
lipticity as usual (e = 1 - b/a, where a and b are the semi-major 
and semi-minor axis lengths, respectively) we found that most 
(S 60%) of the SOAR field objects are not elongated (e S 0.2). 
The counterparts of the INS candidates are also rather point- 
like, with mean ellipticity ~ 0.08 - the most elongated one being 
source XMM J2140 and the two objects possibly associated with 
XMM J0106 (see values in Table0. 



3.2.5. Optical upper limits for XMM J1 046 

The only INS candidate of our sample that remains without an 
optical candidate is source XMM J 1046. The short ESO-VLT ex- 
posure obtained in 2007 revealed no counterparts brighter than 
R ~ 25. The location of XMM J 1046 towards the Carina star 
forming region explains the presence of many bright stars and 
of relatively intense diffuse emission in its optical field. The R 
filter, which includes the strong Ha line, is thus particularly af- 
fected by background intensity variations, impacting the quality 
of the photometry. In the light of this, we asked for additional 
deep SOAR exposures in the V band beginning of 2008. The V 
band excludes some of the strongest lines found in H n regions 
(Balmer emission as well as [N n] and [S n] emission lines) and is 
thus best suited to get a deep upper limit on the optical emission 
of XMM J1046. However, the V filter is still somewhat contam- 
inated by the nebular [O m] emission. 
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Table 6. X-ray spectral analysis of the INS candidates (excluding XMM J 1046) 



Candidate 


OBSID 


t a 


d.o.f 


kT 
lev; 


r \, 

/X,14_ 

{ -1 -2\ 

(ergs cm ) 


c 


Goodness^ 

\ i°) 


r 


/x,14 

(erg s cm ) 


C 


Goodness 

\ '"I 


XMMJ1210 


204710101 


2.1 


22 


^-"-40 


J - -0.7 


14 


4.6 


2.0(4) 


^■"-0.9 


17 


9.9 


XMM JO 106 


150870201 


1.7 


5 


947+116 


2 6 +a9 


5 


18.4 


1.5(8) 


4 6+1.4 

t.U j 3 


5 


19.8 


XMM J0435 


307001301 


7.3 


24 


200^ 


2 00 + " 3 

^•""-0.29 


22 


26.8 


2.4(4) 


2.8(4) 


21 


22.7 


XMM J0314 


201750901 


11.6 


17 


22Tf 
144 + ;ij 
176 +21 


1.7(3) 


24 


76.7 


2.1(4) 


2.5(5) 


21 


60.0 


XMM J2140 
XMM J1259 


008830101 
203020101 


10.5 
9.8 


15 

26 


2 00 +0 - 22 

9 i c+0.26 


17 

44 


58.5 
96.7 


9 Q4+0.24 
9 4^+0.23 


2 59 +0 - 3 

-0.25 

3.1(3) 


9 
35 


6.8 
79.4 


XMM J1250 


303561001 


7.8 


11 


186 + | 


i.ju_ 024 


6 


6.7 


2.6(4) 


1.7(3) 


6 


7.0 



" Exposure times are filtered for background flares. 

* Observed flux refers to range 0. 15 - 3 keV in units of 10~ 14 erg s _1 cm" 2 . 

c The "goodness-of-fit" is derived from a number of 1000 Monte Carlo simulated spectra. 



To test the detection limit specifically in the region covered 
by the error circle of the X-ray source, smalQ sections of the B, 

V and R images were analysed so as to minimize the contribu- 
tions of nearby bright stars and of nebular emission. The frames 
were binned by a factor of 2 in order to increase the signal- 
to-noise ratio. Images of synthetic stars (created using the PSF 
model derived from the data) of progressively fainter magnitudes 
were added at the position of the X-ray source. The composite 
images were then subjected to automatic detection and to mag- 
nitude measurement using DAOPHOT. In this process, the syn- 
thetic star was no longer detected or was rejected (while trying to 
fit the PSF model) when the signal-to-noise ratio was worse than 
S IN ~ 2.5. We then defined the limiting magnitude as the mag- 
nitude of the faintest synthetic star still succesfully measured at 
this confidence level. We give in Table [5] the magnitudes of the 
faintest simulated star (m syn ) detected at a signal-to-noise ratio 
(Tdet, along with the measured magnitude (/«d e t). The discrepancy 
between the simulated and detected values of the R magnitude 
is due to a somewhat poorly defined PSF model in the wave- 
length range most affected by background intensity variations. 
Taking the limiting magnitudes as upper limits on the brightness 
of the optical counterpart of XMM J 1046, Table [5] also lists the 
implied X-ray-to-optical flux ratios corrected for photoelectric 
absorption and interstellar extinction. In order to compute the 
unabsorbed fluxes we adopted the best blackbody fit parameters 
of the source (see Sectio n [4]) and derived Ay = 1.96 using the 
iPredeh l & Schmitt ( ; 19951) relation between the X-ray absorp- 
tion and optical extinction. F or the other bands, w e adopted the 
Ax/Ay extinction relations of lCardelli et alj dl989l) . 

4. X-ray analysis 

XMM J 1046 is a very likely new thermally emitting INS with 
properties that, at first glance, seem similar to those of the M7. 
A detailed X-ray and optical analy sis of this object has recently 
been reported in iPires et all d2009l) . We recall below the main 
results presented in this paper and refer to it for details on the 
reduction and analysis of the X-ray data of this and the other 
seven INS candidates. 

4.1. XMM J1 046 

XMM J 1046 was detected altogether in 16 XMM-Newton obser- 
vations of the well studied binary system Eta Carinae and of the 
Wolf-Rayet star HD 93162 (WR25), as well as in one Chandra 

12 54" x 46" , 128" x 100" and 70" x 70" respectively for the R, B and 

V images. 



observation of the Carina Nebula. Although the observing condi- 
tions were far from optimal - in many observations XMM J 1046 
was located near the edge or in CCD gaps and often the effective 
exposure times were short - we have analysed the available X- 
ray archival data on this source, which span more than six years. 

The event files were reprocessed using SAS 7.1.0 and 
CIAO 4 0.l|H applyin g standard procedure. Using XSPEC 12.4 
(lArnaudlll996l 120041) . we tested different models (blackbody, 
power-law, bremsstrahlung, Raymond-Smith, . . . ). All param- 
eters were allowed to vary freel y. Due to th e low number of 
counts, we applied the C-statistic (Cashl fT979h in order to derive 
the best fit parameters and their uncertainties. Whenever possi- 
ble for a given observation, data from all EPIC cameras were 
analysed simultaneously to better constrain the spectral parame- 
ters. 

We have shown that the X-ray emission of XMM J 1046 is 
best described by a soft blackbody, with mean kT ~ 117 + 14 eV 
and N H ~ (3.5 + l.l)x 10 21 cirT 2 and a stable 0.15-3keV ob- 
served flux of fx ~ (1.03 + 0.06) x 10~ 13 erg s _1 cirT 2 (errors are 
3 <t). There is no evidence for a hard non-thermal component. 
In this work we adopt the notation Fx to denote the unabsorbed 
fluxes while fx is used for the observed flux. The Nh is a fac- 
tor of ~ 10 higher than those typical of the M7 and is consistent 
with the one towards Eta Carinae, which is located at ~2.3kpc 
(ISmithll2006l) . The distance to XMM J1046 is likely to be com- 
parable. Together with the optical data, the high lower limit on 
the X-ray-to-optical flux ratio, \og(Fx/Fy) ^ 3.1*?'?, piratically 
excludes any other possibility than an INS. Timing analysis ap- 
plying a Rayleigh Z 2 test shows no pulsations to a rather non- 
constraining 30% upper limit (3 cr), in the 0.073- 100 s period 
range. 

4.2. Other INS candidates 

Each one of the other INS candidates was serendipitously ob- 
served only once by the XMM-Newton detectors. The event 
files were reprocessed using SAS 8.0.0 applying standard pro- 
cedure as for source XMM J1046. For the spectral analysis, 
the low signal-to-noise ratio spectra were only fitted assuming 
an absorbed blackbody or power-law using XSPEC 12.4. Being 
roughly one order of magnitude fainter than XMM J 1046, the 
low number of counts prevents well constrained spectral fits; in 
particular, it does not permit a strong constraint on the value of 
the column density. Blackbody fits show rather high tempera- 
tures, typically kT <: 200 eV, and column densities ranging from 
to 4.2 x 10 21 cm~ 2 , at 68% confidence level. Although appar- 

13 http : // cxc . harvard . edu/ciao4 . 0/ index . html 
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ently in disagreement with our selection criteria (kT < 200 eV), 
the significant errors in HR explain the inclusion of these X-ray 
sources among the selected sources. We thus decided t o hold A^h 
fixed at the Galactic value dDickev & Lockmanlll99Q, Table [2) 
in order to obtain better constrained spectral fits. The results of 
the blackbody and power-law fits for these sources can be seen 
in Table [6] (errors are 1 cr). The quality of each fit ("goodness") 
corresponds to the fraction of simulations yielding a better fit 
statistic than the actual data, with high values implying bad fits. 

Whereas the relatively large number of counts collected for 
XMM J 1046 safely allows one to exclude a power-law shape for 
its X-ray spectrum, the much lower signal-to-noise spectra of the 
seven sources studied here are equally well described by either 
a hot blackbody or by a soft power-law energy distribution, with 
spectral indexes usually larger than 2. 

5. Discussion 

At energies below 2keV, the X-ray spectra of quasars and 
Seyfert I galaxies usually show an excess relative to the extrapo- 
lation of the high energy (2- lOkeV) power-law. It is generally 
believed that this "soft excess" is produced by the scattering of 
thermal optical/UV photons from the accretion disk surround- 
ing the central black hole to soft X -ray energies, by a po pula- 
tion of ambient hot electrons (e.g. Atlee & Mathur 2009, and 
references therein). Over a limited energy range - in particular, 
that covered by the sensitivity of the ROSAT instruments, 0.1- 
2.4 keV - the energy distribution of the soft exces s is roughly 
descr i bed by a power-law of spectral inde x ^ 2 (e.g. lBoller et alj 
119961; iGrupe et aljfl998l iPageet al.|[2004b . This is in agreement 
with the results obtained in Section |4] for the INS candidates 
with optical counterparts (Table [6). Moreover, X-ray emitting 
quasars cluster in a region of the HRi x HR2 diagram which is 
right above selection region III (contours in Fig. [TJ, the location 
where it is expected to find sources with emission compatible 
with the hottest blackbodies that entered our selection. At faint 
fluxes, the 2XMMp sources show rather large HR errors, reflect- 
ing the fact that the low number of detected photons prevents a 
precise determination of the source spectral energy distribution. 
As it is evident from Fig.Q] several of the INS candidates share 
hardness ratios common to both (galactic and extragalactic) pop- 
ulations. For these reasons, our search procedure is sensitive to 
optically faint AGN, especially those showing soft X-ray excess. 

Magnetic CVs exhibit a low temperature component de- 
tectable in soft X-rays or at UV wavelengths, depending 
whether the source is in a high or low a ccretion state (e.g. 
Ramsav & Cropper 2004; Ramsav et al. 2004). This emission is 
due to the reprocessing of hard X-ray photons (of several tens 
of keV) produced in a shock of the accretion flow into the pho- 
tosphere of the white dwarf. In principle, faint polars showing a 
particularly large soft X-ray excess and with uncatalogued opti- 
cal counterparts can be present in our sample of INS candidates 
as well. 

Active coronae of late-type stars manifest themselves as soft 
X-ray emission with luminosities generally below 10 31 ergs -1 . 
This popu lation dominates soft X-ray samples at low Galactic 
latitudes dMotch et al.ll 19971) . It is well known, however, that the 
high latitude X-ray sky at flux level 10~ 14 ergs' 1 cm" 2 is domi- 
nated by AGN (e.g . IMaccacaro et~aT]|1982t iBarcons et al]|2007t 
MateosetaL|[2008) ; the stellar content, although more important 
in soft X-ray-selected samples as ours, should be less than 10%. 
In the optical, none of the INS candidates show colours typical 
of M, G, K stars (Fig. [5j. In any case, it is expected that late- 
type stars have log(/ x //R) < -1 (e.g. IMaccacaro et aT] fl982; 



IStocke et al.ll 199 lb which is clearly not the case for any of our 
INS candidates, all with log(/x//R) > 0. Therefore, we do not 
expect any of the investigated X-ray sources to be identified with 
late-type stars. 

Based on the X-ray/optical associations discussed in 
Section 13.2.31 the objects found inside the error circles of the 
INS candidates are very likely to constitute their true optical 
counterparts, all with a very low (%5%) probability of chance 
association, with the exception of XMM J0106. The strong UV 
excess exhibited by several of the sources (XMM J 1259, XMM 
J2140, and XMM J1250) is clearly consistent with those of AGN 
and CVs (Pig.|5j- It is worth noting, however, that optical colours 
alone would not be sufficient to distinguish these X-ray sources 
from the population of (X-ray and radio-quiet) blue objects lo- 
cated at high Galactic latitudes. In spite of the more average 
colours of the optical candidates for sources XMM J1210, XMM 
J0435 and XMM J0314, their high X-ray-to-optical flux ratios 
leave no doub t that they are likel y of extragalactic origin. Still 
according to lBarcons et al.1 (120071) . the fraction of X-ray sources 
with log(/x//R) > 1, where obscured AGN are expected, al- 
though more important in hard X-ray-selected samples, is of just 
a few percent in soft ones. As mentioned in Section 13.2.31 two 
out of seven of our candidates show log(/x//R) > 1, XMM 
J0106 and XMM J2140. These are the sources with the most 
elongated optical candidates as well. In particular, the two opti- 
cal objects found in the error circle of source XMM JO 106 cor- 
respond to the most dubious of the X-ray/optical associations. 
They have been detected only in the R filter, preventing colour 
estimation, whereas the optical candidate of source XMM J2140 
has been well detected in the U, B and R filters. XMM JO 106 
was detected in a very short exposure and therefore its HR errors 
are worse determined; however, the two possible optical counter- 
parts define an already very high log(/x//R) ~ 2, characteristic 
of the most extreme clas ses of X-ray emitters as AM Her sys- 
tems and BL Lac objects ( Schwo pe et al.|[T999l) . and it is worthy 
of further investigation. 

Our best INS candidate is XMM J 1046, for which no optical 
counterpart is detected down to the limiting magnit ude of our 
presen t data. On the basis of the analysis reported in Pire s et ail 
(120091) . the identification of XMM J 1046 with an INS is very 
likely. This could be the first example of a presumably radio- 
quiet and X-ray dim INS, located at a significantly greater dis- 
tance than the M7. Its spatial location in the Galactic plane and 
the derived value of the column density suggest that it may be 
physically associated with the Carina Nebula, a giant H 11 region 
harbouring a large number o f massive stars and wit h ongoing 
active star formation (see e.g. ISmifh & Brooksll2007l and refer- 
ences therein). 

Several attempts aiming at identifying new thermally emit- 
ting INSs have been carried out in the years that followed the dis- 
covery of the seven sources, usually by cross-cor relating ROSAT 
data with optical, radio and IR catalogues (e.g. Rutledge et al. 



2003 
120071) . 



lAgiieros et alj|2006l IChieregato et al.ll2005t iTreves et al 
However, these searches are hampered by the large 
ROSAT positional errors at faint fluxes, especially in the pop- 
ulated regions of the Galactic plane. For these regions, many 
spurious (low significance) candidate optical/IR counterparts en- 
ter the X-ray error circle of a given ROSAT source. This makes 
the probability of erroneously assigning an identification signif- 
icant, thus excluding the source as a potential INS candidate. 
The final effect is that the Galactic plane is largely "avoided" by 
usual cross-correlation algorithms (see e.g. Rutledge et al. 2003, 
Fig. 1), and proposed INS candidates are usually located at high 
Galactic latitudes. 
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A long-term project which investigates INS ca ndidates se- 
lected from the ROSAT Bright Source Catalogue (Voge s~et~aT] 
119991) . making use of follow-up investigations with the Swift 
satellite, is presently bei ng conducted by R. E. Rutledge, D. 
B. Fo x and collaborators (Rutled ge et alj|2008t iLetcavage et alj 
120091) . Follow-up X-ray and optical observations of one candi- 
date with no evident counterpart led to the discovery of Calvera, 
a likely INS (which exact nat ure is, however, still u nclear) with 
a very large F X /F Y Z 8700 dRutledge et al.ll2008l) . Not an ex- 
ception, Calvera is located at high b, implying a remote distance 
of d = 8.4 kpc and a vertical velocity in excess of ~ 5 100 km s _I 
to explain its current position well above the plane, considering 
standard cooling time and thermal emission similar to the M7. 
Instead, it is more likely that Calvera is a nearby radio pulsar. 

According to the de tailed population synthesis calculations 
of iPosselt et alj ((2008), it is expected that, in general, new 
unidentified cooling INSs at faint soft X-ray fluxes are to be 
found in the Galactic plane, at small angular distances from 
their birth star forming regions - in particular, of the rich OB 
associations located beyond the Gould Belt, as Carina, Vela, and 
Cygnus-Cepheus. At the flux limit applied in our search, the ex- 
pected number of cooling neutron stars having the same prop- 
erties as those exhibited by XMM J1046 - i.e. a slightly hot- 
ter source located at a greater distance relative to the M7 - is 
of only ~50-80 in the whole skjQ, which translates into less 
than one source to be present in the 2XMMp catalogue. The fact 
that we found one source exhib i ting th ese unique characteristics 
only corroborates IPosselt et all (|2008) conclusions that a search 
for new cooling INSs should not be blind, i.e. one has to look 
for preferentially in the most promising regions of the sky - i.e. 
near more remote OB associations. 

Population synthesis constitutes a fundamental tool in order 
to compare observational results from e.g. surveys with theo- 
retical expectations for a given population of astrophysical ob- 
jects. If a large number of objects are known, as for the case 
of radio pulsars, it allows testing many of the population phys- 
ical (and theor etically unknown) properties. Re cent examples 
are the work of iFaucher-Giguere & Kaspl d2006l) . conducted for 
the population of radio pulsars, and the study carried out by 
Stor y et alj d2007l) with the population of disk millisecond radio 
and y-ray pulsars. 

In the particular case of thermally emitting INSs, population 
syntheses conducted in the early 90's predicted that a large num- 
ber of old neutron stars accreting from the ISM would be de- 
tected bv ROSAT dTreves & Colpj[l99Tt iBlaes & Madaull 19931; 



lMadau&Blaeslll994l) . It soon became clear that the unrealis- 
tic overprediction was due to several simplfying assumptions 
used in the model. In fact, several mechanisms act in order to 
inhibit accretion such as the presence of the magnet ic field of 
the neutron star (Toropinaet al. 2003; Ikhsanov 200 7]) and that 
of an inf ailing material which is weakly magn etized dPerna et alj 
120031) or heated by the emergent X-rays (e.g. Blaes et al. 1995). 
In any case, it is believed that accretion should proceed at a rate 
well below the Bondi-Hoyle one, thus affecting the actual num- 
ber of sources that can be observed by current X-ray missions. 
The population of cooling ne arby INSs from the Gould Belt has 
been extensively modeled bv lPopov et a l. (2000, 2003, 2005D. In 



14 The rang e shown here for th e expected number of sources is 
derived from Posse lTet alj J2008T) adopting the dotted and hashed 
curves of their Fi g. 5 as lower and upper limits, respectively. See 
IPosselt et al] l200§ for details on the assumed parameters. EPIC pn 
counts were converted to ROSAT PSPC counts using WebPIMMS 
(http : //heasarc . gsfc .nasa . gov/Tool s/w3pimms . html). 



iPopov et"a"fl (120061) . the expected number of sourc es at a given 
flux w as used to constrain cooling curves of INSs. IPosselt et all 
(2008) improved many of the features of their population syn- 
thesis model, mainly by providing a new description of the ISM 
distribution and that of the massive progenitors. Interestingly, 
they used the spatial distribution of the synthetic population in 
order to define search strategies to identify new cooling INSs. 

The expected number of such sources to be detected by cur- 
rent X-ray missions is strongly model-dependent. It involves 
from the parametrization of properties inherited at birth (kick 
velocity and mass distribution of the compact remnant) to that 
of the cooling history and emission mechanism of the neutron 
star. Together with realistic descriptions of the Galactic gravita- 
tional potential, the ISM distribution and the detector character- 
istics, these ingredients determine the neutron star thermal and 
dynamic evolution and detectability. Similarily, from the obser- 
vational point of view, selection biases and the area effectively 
covered by the survey (considering the constraints applied while 
correlating X-ray sources with other wavelengths) have to be 
taken into account in order to derive correct observational up- 
per limits on the number of sources at a given flux. Overall, even 
the simplest models end up having a level of complexity which is 
made evident from the number of assumed hypotheses - many of 
them involving significant theoretical and observational uncer- 
tainties. We are currently developping a model for this Galactic 
population which aims at making the least number of uncertain 
assumptions - even at the expense of having a less sophisticated 
description of the real population - while, based on Monte Carlo 
simulations and leaning on what it is observed, in particular, with 
the XMM-Newton, would be statistically robust in order to con- 
strain important properties of this population of unique Galactic 
remnants especially beyond the Solar vicinity. 



6. Summary and conclusions 

Our search for new thermally emitting INSs in the 2XMMp cat- 
alogue has revealed a handful of interesting and previously un- 
known soft X-ray sources. Deep optical imaging revealed likely 
optical counterparts for six of them which, based on blue optical 
colours and X-ray-to-optical flux ratios around 1-10, identify 
them with, most likely, AGN and CVs. Source XMM J0106 has 
no evident optical candidate and should be further investigated. 

Source XMM J 1046 is a newly discovered thermally emit- 
ting INS, possibly radio-quiet and similar to the M7. As it is ex- 
pected for cooling neutron stars, its X-ray emission is described 
by an intrinsically soft and thermal energy distribution, stable 
on long time scales and with no evidences for magnetospheric 
activity. The present lower limit on the X-ray-to-optical flux ra- 
tio already excludes standard classes of X-ray emitters. Its likely 
location in the Carina Nebula is in agreement with the most up- 
to-date expectations of population synthesis models. 

This work confirms the use of the XMM-Newton catalogue 
of sources as an efficient tool in order to identify new thermally 
emitting INSs and other particular classes of soft X-ray sources. 
A similar search using the Chandra data with its excellent as- 
trometric accuracy would be instrumental to select sources es- 
pecially located in the Galactic plane, where it is more likely to 
find new INS candidates. 
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